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Several cases of wave motion without damping were solved to obtain 
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impulse loading and a square wave stress impulse loading for t ^ aJ c 
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CHAPTER I 


Introduction 

Although there have been numerous analyses published on wave 
motion from many standpoints, no publications on wave motion with 
coulomb damping have been uncovered. This treatise is an initial 
attempt to analyze wave motions with coulomb damping, using Laplace 
transforms as an aid in solving the partial differential equations. 
Wave motion without coulomb damping is presented first to establish 
the fundamental characteristics of wave motion. Then, the solutions 
which were successfully obtained of wave motion with coulomb damping 
are presented. 

in the cnapcer on wave motion witiiuut coulomb damping, semi- 
infinite and finite rods are considered. For the semi- inf inite rod, 
three types of impulse loadings are analyzed: step stress, square 

wave stress, and step velocity. The finite rod is analyzed for the 
right end boundary conditions of free and fixed. Impulse loadings of 
step stress, square wave stress, square wave velocity and triangular 
wave stress are considered for the free right end condition, and 
square wave stress and triangular wave stress for the fixed right end 
boundary condition. 

In the chapter on wave motion with coulomb damping, only a few 
cases are solved because of the difficulty in obtaining solutions; 
however, both the semi- inf inite rod and finite rod are considered. 


For the semi- infinite rod, impulse loadings of step stress and square 
wave stress for t = a/2c and t 5 a/c are analyzed, and for the finite 
rod, a step stress is analyzed. 

In Chapter V the effects of friction- induced wave motion are pre- 
sented for the inducing wave propagation velocities of twice, equal to, 
and one-half of the propagation velocity of the induced wave. In the 
last chapter, the advantages and disadvantages of the Fourier trans- 
form technique as applied to wave motion are presented. 


V 
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CHAPTER II 
Previous Work 

Although almost every book dealing with dynamics, vibration, elas- 
ticity, mechanics or physics devotes some space to wave motion without 
damping, nothing seems to be available on wave motion with coulomb 
damping. Some of the works used as background information for the pre- 
paration of this analysis were Jacobsen and Ayre 1 , Timoshenko 2 , 

Burton 3 , and Tong 4 . The technique of solving the wave equation using 
Laplace transforms was obtained from Wylie 5 and Norwacki 6 . The method 
of variation of parameters for solving some of the nonhomogeneous dif- 
ferential equations was obtained from Kells 7 . The Laplace transforma- 
tions were taken from the tables and principles of Tse, Morse and 
Hinkle 8 and from the tables of the CRC Standard Mathematical Tables 9 . 

Some of the material reviewed as background information on 
Fourier transforms was Sneddon 10 , Hildebrand 11 and Erdelyi, et al 12 . 
Transformation tables of Fourier cosine transforms were obtained from 
Sneddon 10 and Erdelyi, et al. 12 . 

Since no material was available on wave motion with coulomb damp- 
ing, this work was generated from the collective background information 
of the bibliography. However, some work has been done on wave motion 
with various types of internal damping. Narwocki 6 and Goldsmith 13 are 
excellent examples of those who have performed analyses of wave 
motion with visco-elastic damping. 
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CHAPTER III 

<• 

Wave Motion Without Damping 

Wave motion is considered in this report only for a continuous 
elastic system which is covered by partial differential equations. 

The media are assumed to follow Hooke's law and are homogeneous and 
isotropic. The treatise will be limited to one-dimensional wave 
motion in prismatic bars where the length of the wave is large com- 
pared to the cross-sectional width. The cross-sectional planes are 
assumed to remain plane and for this chapter, damping is neglected. 

A. Wave Equation Formulation » 

Consider the prismatic bar in figure 1. Let T be the force 

acting on the cross section at some position x, and at x + dx the 

force is T + |^ dx. Using Newton's law of motion, we obtain 
ox 

I P ext ■ ““ 

-T + T + dx = p A dx , ( 1 ) 


4 


where p is the mass density, A is cross-sectional area, and u is the 
displacement. The force at the cross section is proportional to the 
strain (Hooke 1 s . law) ; i.e.. 


T 



where E is Young 1 s Modulus of Elasticity. For a prismatic bar that is 
homogeneous and isotropic, 


S ~ AE 


and substituting into equation (1) yields 


d 2 u £ 5 2 u n 

’ E 


5 2 u 1 d 2 u _ n 

dx 2 c 5 St 2 * 


( 2 ) 


where c = n/e/ p. The general solution to equation (2) can be expressed 
as 


u = f (ct - x) + g(ct + x). 


and regardless of the functions f and g, the argument (ct ± x) leads 
to the differential equation upon differentiating and substituting. 
Then, for u = f(0), 


c 





which demonstrates that c is the wave propagation velocity. 


B. Semi- Infinite Rod 


This rod begins at x = 0 and extends in the positive x direction 
to infinity. It is an ideal rod used for mathematical purposes to 
give wave motion that is undisturbed by a boundary. 

1. Step Stress Impulse Loading 

The Laplace transform technique is used to obtain a solution 
to equation (2) because the solution for transient type impulse is 
desired, and these solutions are easier to obtain using Laplace trans- 
forms. To solve equation (2), two boundary conditions and two ini- 
tial conditions must be known. For example, let the displacement, 
u (x, t) , and the velocity, u t (x,t), be zero at t = 0. Also, let a 
step stress impulse loading be applied to the semi- inf inite rod. 

This means that the two boundary conditions are as follows: (1) 

u (0,t) = — H(t) where T is the applied load, and H(t) is the 
Heaviside Function (when the argument t > 0, H(t) = 1 and when t < 0, 
H(t) = 0) and (2) as x — » oo, u(x,t) is bounded. Taking the Laplace 
transformation of equation (2) with respect to t yields 

d 2U ^ X ^ - - ^ U(x,s) + u(x, 0) + u (x,0) = 0, 
ox c c c 


where s is the new variable. Since the initial conditions are zero, 
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This equation is an ordinary second order differential equation that 
has a solution of 


s s 

— x - — x 

c c 

U(x,s) = Bje + B 2 e 


Since as x -» °°, u(x,t) is bounded, then Bi = 0; therefore, 


U(x,s) = B 2 e 


4 


5U (x,s) 
Sx 



Applying the other boundary condition, i.e., 




u„(0, t) = H(t) or U (0, s ) = fir , 

A. JA£j Afja 


yields 


T 

AEs 



or 



U < x ’ s > = AEs^ 6 

Taking the inverse transform gives a solution to the partial differ- 
ential equation. 


u(x,t) -g (t -f) H(t -f), 


* 
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Essential to the understanding of wave motion is the velocity and 
strain distribution of the rod. These can be obtained by taking the 
partial derivatives with respect to t and x; i.e., 


i 


5u 

at 


Tc , x N 

= 7^ H(t - -) 
AE v c 


bu 

bx 



The distribution of displacement, velocity and strain for various 
times , t i5 is shovm in figure 2. A step stress input yields a step 
velocity and strain but a ramp-type displacement. Also, the dis- 
placement cannot be discontinuous like the velocity and strain, 
because a discontinuity would indicate a break in the rod. 

2. Square Wave Stress Impulse Loading 

Let the initial conditions be zero. At t = 0, a stress is 
applied at x = 0 until t « t. Then the stress is zero at x = 0; 

1 « 6 • j 


> 


U (0,t) = - ^ [H(t) - H(t - t)]. 


( 3 ) 


Equation (2), the differential equation to be solved, has the general 
solution of 


s s 

- x - - x 

U(x,s) = Bxe C + B 2 e 


For a semi-infinite rod, Bi = 0, and 


s x 

u x (x,s) = c e C (4) 

Taking the Laplace transform of equation (3) and substituting into 
(4) for x = 0 yields 


T 

AEs 


(1 




-TS 

e 


) 


U(x,s) 



(1 - e' TS ) 


£ 

c 


X 


u(x,t) 


Tc 

AE 


(t - |) H(t - |) 


(t - T - |) H(t - T - |) 


rin 


Tr 


st = tz L H(t - p - H ( t - x - 


Su = _ _r_ 

Sx ' AE 


H(t - 2) - H(t - t 



The distribution of displacement, velocity and strain is shown in 
figure 3. The square wave stress input yields a square wave velocity 
and strain but a ramp wave front on the displacement tohich is constant 
after t = r. 

3. Step Velocity Impulse Loading 

For a semi- inf inite rod with a step velocity at x = 0, the 
end condition is 

u t (0,t) = vH(t) . 


(5) 
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With the initial conditions zero, the solution to equation (2) is 


U (x, s) = B 2 e 


s 

- — x 

c 


( 6 ) 


Taking the Laplace transform of equation (5) with the initial condi- 
tions equal to zero and substituting into equation ( 6 ), for x - 0 , 
we obtain 


^ = B 2 


s 

- — X 

, .. v c 

U(x,s) = 72 e 


u(x,t) = v (t - H(t - 
u t (x,t) = vH(t - f) 
u x (x,t) = " c H ^ t " c^ • 


These distributions are shown in figure 4. It is interesting to 
note the simularity of the distributions for a step stress and a 
step velocity. 

This section illustrates the characteristics of wave motion 
in its simplest form to form a basis for the analysis of the more 
complex wave motion which is to follow. An understanding of these 
simpler forms simplifies the understanding of the more complex and 
allows the individual to estimate the behavior of the more complex 


wave motion. 


% 
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C. Finite Rod 

Although the finite rod has a few of the characteristics of the 
semi- inf inite rod, the finite rod exhibits unusual behavior at the 
boundaries caused by the added phenomenon of reflection. 

1 . Right End Boundary Condition Free 

The term "free" here means that the right end is unsupported, 
and the stress at the end is zero unless a stress is specifically 
applied. 

a. Step Stress Impulse Loading to the Left End 

A rod of length, £ 9 has a compressive step stress applied 
when t = 0 and at x = 0; therefore, the boundary conditions are as 
follows : 


u (°» t > = - TZ net) or u (O.s) 

X A. Hi X 


AEs 


( 7 ) 


u (i,t) = 0 or U (£,s) = 0. 

X X 


( 8 ) 


The general solution to equation (2) for zero initial condition is 


¥ 


¥ 


U(x,s) 


s 

— X 

= B^ 0 + B 2 e 


s 

— x 

c 


; 


U (x,s) 

X 




Substituting equation (7) with x - 0 and equation (8) with x 
gives 


Tc 

AEs 


2 ” Bi - B 2 


0 = e B x - e B 2 


n 


Bi - 


Tc e 


AEs 2 1-e 


U 


and B 2 = 


Tc 


AEs 2 1-e 


2& 


U(x,s) = 


Tc 


AEs 2 1-e 


2J& 


-| s (- ^ + |)s 
i + e 


Using the binomial expansion: yields 


U c -Ms - ^ 


2JI 


= 1 + e 


+ e 


+ e 


+ . . 


1-e 


Therefore, 


U(x,s) 



+ e 


(.21 . £ )s (- il 

v c c J v c 

+ e 



(- 

+ e 


6£ 

c 


(continued on next page) 
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r- (- M. + *) s (- + i)s (- + £)s (- ^ + -)s 

1 v c <r ( v c c' , v c c' , v c c , 
e + e +e + e + 


and taking the inverse Laplace transform yields 


u(x,t) 








+ 



.ii + sN + ft-ii + a 

c c J \ c 



The velocity equation is 


u t (x,t) 



- ) + H 





+ H 


t _ M _ * 


H 


- ^ + S] + H 

c c. 


t - 


4,g 




(9) 



+ H 


( 10 ) 
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and the strain equation is 


u x (x>t) = a* 


H 


- - ) + H 
c 


- — + H 


c 


+ H 



H 



+ H 





(ID 


ghown in figure 5 is a representation of the displacement, strain and 
velocity distributions at time, t = 0. At t = 0, the contributions of 
each term of the above equations have not yet reached the rod; however, 
as time increases, the displacement, velocity and strain waves move in 
the directions of the arrows, and some of the waves immediately enter 
the domain of the rod. Figure 5 is used as a step to obtain figure 6. 
However, similar figures to figure 5 will not be shown in the solution 
to other types of examples to follow, but should be understood as a 

step that was taken but not shown. 

To simplify the solutions which will follow, the equations will 
be written in summation notation. Equations (9), (10) and (11) in 
this notation would be as follows: 


00 




H t - 
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■ i {i[* ^ - 1)] + z [»g - ^ + 1)] } • 

n=0 n=l 


“*<*•'> ' • a {Z [ H (‘ ' - f)] - 1 [ H (' - ¥ + !)] } 

n=0 n=l 


It is advantageous to make some observations about fig- 


i i 


are 6 at this point. Notice that the displacement ana velocity con- 
tinue to increase while the strain fluctuates. The velocity at first 
increases from the left where the load is applied due to a compressive 
wave, then the velocity doubles at the right end and increases from 
the right due to an expansive wave. The strain reaches a constant 
value throughout the rod when the displacement is of a ramp shape 
over the length of the rod. The strain is zero over the length of the 
rod when the displacement is constant over the length of the rod. 

b. Square Wave Stress Impulse Loading to the Left End 

The square wave input is more realistic in practice than 
the step input, because at some time^ x, after the stress has been 
applied, the stress input ceases; therefore, the rod with a square 
wave input does not continue to increase its velocity as it did with 
a step input. Considering the same rod as before with zero initial 
conditions and 


^(O.t) ' ae 


H(t) - H(t - T ) 


or V°’ s) = aeF <e " TS • 


u (Z, t) = 0 or U (i,s) = 0, 

X X 


the solution to the equation of motion (2) is 


00 



n=0 




with the corresponding velocity and strain distributions 



( 12 ) 



The velocity and strain distributions are shown in figures 7(a) and 
(b) for t = The velocity wave travels down the rod and doubles 
at the right end. As the wave returns, it reduces to its original 
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value. When it reaches the left end, the same wave motion occurs as 

* did at the right end. The strain wave travels along the rod until it 
reaches the end, and then both ends of the wave go toward the middle 

' of the wave until the wave vanishes at ct/4 from the end of the rod. 

After the wave vanishes, it reappears at the same position with the 

opposite sign, and the wave expands until it has a width of t. When 

it is the width of x, the strain wave has reached the end of the rod, 

and the wave moves in the opposite direction from which it initially 

moved. When the wave reaches the opposite end (x = 0), the same wave 

motion occurs at that end as did at the other end (x = Z) . 

Figure 8 illustrates the displacement distribution of the 

0 

square wave input of pulse duration, x = ■^. When the wave reaches 
the end of the rod (x = Z) , the wave form changes. This is caused by 

A 

the interaction of the two ramp-shaped waves meeting at the end of the 
rod. To give a further illustration of the displacement of the rods, 
the displacement versus time is shown in figure 9 for three rod posi- 
tions x « 0, Ji/2 and Ji. A rod when struck by a square wave of pulse 
Ji 

duration, t = — does not move en masse, but only portions of the rod 
move at one time, like a worm crawling. The ends of the rod move in 
a similar manner with the end, x = i, at a phase lag of t = £/ c to the 
end, x = 0, while the middle moves only 1/2 the distance at one time, 

* but the times occur twice as often. The reason the ends move twice as 
much as the middle is because a compressive wave moves down the rod tc 

^ the end, x = Ji, then converts to an expansive wave of the same magni- 

tude, thereby causing a displacement twice that due to the original 



wave. This wave motion is not typical of nature in that nature seldom 
applies a stress of wavelength t = This will be illustrated next 

with a square .wave velocity impulse loading. 

c. Square Wave Velocity Impulse Loading to the Left End 

This example requires some forethought before the desired 
solution can be obtained. If we wish to simulate the striking of a 
golf ball, a baseball, or some other similar action, the duration of 
the square wave is important. This pulse duration is generally the 
time that it takes the pulse to go the length of the elastic body and 
return; i.e., t = — . Other pulse durations may be specifically 
applied, but in general, they are not true simulations of nature. To 
illustrate this point, a solution will be obtained with a general 
pulse duration width of t. The difficulty at this point is defining 
the boundary condition at the end where the square wave velocity is 
applied. When the velocity is applied, the end condition is defined 
as 


u t (0,t) = vH(t) or U(0,s) = , 


where v is the applied velocity, and when the velocity is released. 


the end condition is 


V 
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P 


It is difficult to apply this end condition in this form; however, if 
the end condition could be totally defined in terms of u x , then it 
could be easily applied. This can be done by obtaining the solution 
of a rod which is free on the right end and has a step velocity impulse 
loading on the left end; i.e., let the initial conditions be zero, and 
the boundary conditions be 

u fc (0,t) = vH(t) or U(0,s) = jg 

u x U,t) *0 or u x U,s) = 0 . 

The solution to the equation of motion, equation (2), is 


n=0 


s* + f 


n=l 


The corresponding velocity and strain equations are 




n=0 


n=l 


§ 




n=0 


(14) 


n=l 


These are shown in figure 10. The velocity, v, travels down the rod 
and doubles at x * I. As the rod attains a velocity of 2v, the strain 
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is being reduced to zero. When the wave is returned to the origin, 
the strain changes sign and the velocity is reduced back to v. This 
is because the boundary conditions require the velocity to be v at 
x = 0. This is as though the input is attached to the rod. If the 
input had not been attached to the rod, then the rod would have left 
the input at a velocity twice the input velocity. The fact that the 
strain changed signs indicated that the input was no longer pushing 
but was pulling. 

The objective of this exercise was to obtain an end con 
dition at x = 0 entirely in terms of u x . Therefore, using equation 
(14) yields 


u(0,t) 

X 


v 

c 


oo 00 

y (-D n h u - t 1 ) + X ( " i)n h G 



n=0 n=l 

which is the boundary condition at x = 0. To limit this equation for 
t g 2 jg/c, we eliminate the summation and the latter Heaviside func- 
tion. Thus we obtain the left end boundary condition while the 
velocity input is applied: 


u x (0,t) 


v 

c 


H (t ) . 


To obtain the total left end boundary condition, this equation is 
combined with u x (0,t) = 0 (after the velocity input is removed) to 

ob tain 

u x (0 ,t) = - J jjH(t) - H(t - t) 


for t i . 
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By using this boundary condition with the right end boundary condi- 
tion, 

U x u>t) = 0, 

and setting the initial conditions equal to zero, the solution to the 
equation of motion (2) is 


00 



for x ^ 2$/c. For t = 2^/c which is generally characteristic of 
nature, the equation reduces to 



and the corresponding velocity and strain equations are 



Figure 11 shows the velocity and strain distributions. The velocity, 
v, and strain, -v/c, traverse the rod, simultaneously. When the wave 
reaches x = i, then it is reflected with the velocity equal to 2v, 
while the strain is reduced to zero. When the wave reaches x = 0, 
then the rod leaves the input at a velocity of 2v and all wave motion 
stops. To explain this phenomenon in terms of a golf club striking a 
golf ball, we would need to assume that the club was rigid and did 
not slow down while striking the ball. When the club strikes the 
ball, a wave goes across the ball bringing the velocity of the ball 
to the velocity of the club and putting a strain in the ball. As the 
wave is reflected, the strain is relieving itself by pushing against 
the club, thereby causing the velocity of the ball to reach twice 
that of the club. It is now apparent why x = 2&/c is characteristic 

of nature. 

To get further information concerning the step velocity, 
let x = i/2c, then equation (15) becomes 



(4n-H)l x\ 
2c cj 


• H 


(4n+l)i? 

2c 




H (t 



e - + f) = (t - | 



with velocity and strain equations of 



H t 


(4n+l)l x\ 
2c c J 


00 
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n=l 



H t 


(4n+l)l 

2c 



00 



n=0 



H t 


(4n+l ) Z x 
2c " c 




(4n+l)l 

2c 



These equations are the same as that of the square wave strain input, 
equations (12) and (13), with the exception of the constants in front 
of the equations. Therefore, figures 7(a) and (b) are illustrations 
of the velocity and strain distributions, and only the scales need to 
be changed. Figure 8 is also applicable to the displacement. 

d. Triangular Wave Impulse Loading to the Left End 

Consider the triangular wave shape shown in figure 12. 
The wave shape may be expressed by the equation 


f(t) = T ^jtH(t) ' * * &)+(' -■£) « (t - £)} . 


Therefore, the boundary conditions are 


“*<"•«> - - - 2 (' - £) H ( c - £) + (‘ - i) H ( c • £ 


u x (i,t) = 0 . 

With the initial conditions set to zero, then the solution to the 
equation of motion is 


u = 


2Tc" 

AE$ 


t - 


2nl _ x\ f t _ 2nf _ x 

c c ) V c c 


* - 2 t 


( 8 n+lH 

4c 


n=0 


H 




( 4 n+l)l 

2c 


► I [(* - ♦ i) H (‘ - ¥ - !) - 2 (« - ^ + i) 

n=l L 

- + f) - (« - - I/ H (t - + 


H t 


with a velocity and strain of 


_ 4Tc‘ 

U t " AEf 




c c 


c c 


.tU.sagii 


n=0 


H t 




(equation continued on next page) 


° olx ^ lx, k° I* 
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00 



„ _ (8n+l)l ± x'\ 1 /' k (4n+l ) j, , ^ TT (4n+l)i , *V 1 

“ ^ 4c ’ cj ' V" " 2T 2c T c/J )' 


The only difference between the velocity and strain equations is the 
factor, c, and a sign change in the first group of terms. Figures 13 
and 14 reveal the effects of this difference. The interesting thing 
about these figures is the characteristics of the wave motion at the 
boundaries. The velocity doubles at the same time the strain changes 
sign. The displacement is shown in figure 15, and continues to 
increase as expected. Notice the unusual wave front and its behavior 


at the boundaries. 
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2. Right End Boundary Condition Fixed 

The fixed right end boundary condition implies that the end 
is rigidly clamped. The left fend is allowed to move while the fixfed 
end can have no motion. Thus, the impulse loads will be applied to 
the left end only. 

a. Square Wave Stress Impulse Loading to the Left End 
Of interest in this example is the wave shape at the 
fixed end, x = &. In solving this example, the initial conditions 
will be set to zero. The boundary conditions are 


u x (0,t) - - jg [ (H(t) - H(t - t)] or U x (0,s) 



u(j£,t) =0 or UU,s) = °- 

Using these boundary conditions, the solution to the equation of 
motion (2) is 
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with the corresponding velocity and strain being 


00 



= x /v t.ii a r.»4 


i/-' 

n=0 


2*1 . j. u (+ 

c c J ' “ V 


2nl _ xM 
c ' c/J 


OO 



The velocity and strain distributions are shown in figures 16 and 17 
for t = Hi 2c. At the fixed end, the velocity changes direction, and 
the strain doubles itself; and at the free end, the strain changes 
direction, and the velocity doubles itself. The displacement is 
shown in figure 18. Notice the double slope when the wave reaches 
fixed end. 

b. Triangular Wave Stress Impulse Loading to Left End 

Consider the triangular wave shape shown in figure 19. 
This wave shape may be expressed by the equation 


f(t) = T (1 / t) [ (t - t) H(t) - (t - t) H(t - t)]. 


Therefore, the boundary conditions are 


u (0,t) = - T” [(T - t) H(t) - (t - t) H(t - t)] 

X Aij T 

u(jg,t) = 0. 

With the initial conditions set to zero, then the solution to the 
equation of motion (2) is 



with the corresponding velocity and strain 



n=l 


(equation continued on next page) 
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The velocity and strain distributions for t = jg/ 2c are shown in fig- 
ure 20. Notice at the fixed end the way that the velocity changes 
directions, and the strain doubles itself. At the free end the 
exact opposite occurs; i.e., the strain changes direction in the 
same way the velocity did at the fixed end, and the velocity doubles 
in the same manner as did the strain at the fixed end. The 
displacement is shown in figure 21, 


CHAPTER IV 


Wave Motion with Coulomb Damping 

The wave motion with coulomb damping follows the same laws as 
wave motion without damping, with the exception that the wave motion 
does not continue indefinitely but is damped out. The difficulty in 
analyzing wave motion with coulomb damping is maintaining the damping 
force in the opposite direction to the velocity and at the same time 
having a linear relationship. One reason for this is that in coulomb 
damping, the damping force is not a function of velocity but is a 
constant; therefore, when the velocity changes direction, no non- 
linear terms can exist in the damping term which will correspondingly 

change direction. 

A. Formulation of Wave Equation with Coulomb Damping 

Consider the prismatic bar of figure 22. Using Newton's Law, 
we obtain 



‘ ext 


= mu. 


and substituting the forces from figure 22 yields 
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where sgn[u] is positive if u is positive and negative if u is nega- 
tive. Sgn[u] is a nonlinear term, and the objective will be to 
define a suitable substitute for sgn[u] that is linear. This will 
not be determined until a particular problem is decided, then depend- 
ing on the direction of the velocity, the appropriate substitute will 
be used. Since 


T = AE 


5u 
dx 5 


then 


AE 



dx - pA dx 



sgn[u] F dx 


d 2 u _ £ d 2 u 
dx" E St 2 


sgn[i] — . 


Since the wave velocity is defined by 


then 


c 


2 



5 2 u 1 d 2 u 

dx 2 7 s St 2 


sgn[u] 


AE ’ 


( 16 ) 


which is the wave equation including coulomb damping. 


B. Step Stress Impulse Loading to a Semi-Infinite Rod 


Like the previous chapter, Laplace transform techniques are 
used to solve the wave equation, and two initial conditions plus two 
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boundary conditions need to be satisfied. The wave equation (16) 
must be altered to fit the particular example (see figure 23); i.e., 
the coulomb force per unit length, F, is applied only over the length 
of the bar that has motion. Therefore, 


5 2 u 1 d 2 u _ _F_ x. 

dx 2 “ c 2 5t 2 AE aK c 


(17) 


Taking Laplace transforms and setting the initial conditions to zero, 
we obtain 


5 2 U (x,s) 



U(x,s) 



AEs 


(18) 


The solution to the homogeneous equation is 


s s 

- x - ~ x 

U c (x,s) = Bje + B 2 e 


To obtain the particular solution, substitute 


s 

— x 
c 

U (x,s) = Dxe 

P 


into the differential equation (18): 


-2 




+ D 


s 2 


xe 


s 

— x 

c 


s 2 

—5 Dxe 




AEs 


D 


Fc 

2AEs y 
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U (x,s) 

H 


Fcx 

2AEs 2 e 


s 

— X 

c 


U(x,s) = U c (x,s) + u p (x,s) 


s 

— X 

c 

B x e + B 2 e 



Fcx 

2AEs 2 6 


S_ 

c 


X 


The boundary conditions for this example are 


u( 0 ,t) 



or U(0,s) 


T 

AEs 


As x oo u(x,t) is bounded or U(x,s) is bounded. Using these 
boundary conditions to solve for the constants B x and B 2 > the solu- 
tion becomes 


U(x,s) 


Fc 2 

“ " 2AEs 3 6 




Fcx 

2AEs^ 6 


s 

“ X 

c 


Taking the inverse transform yields 


u(x,t) = — 


T(t - f) - f (t - f ) 2 


— (t - — 

2 ^ c' 


H(t-f). 


2T 

Let x Q = ct and a = — . Then, 


u = 4 ^ (x - x 0 )(x° + x - 2a) H(t - |), 


(19) 


with the corresponding velocity and strain being 


“t ■ (a - X o) H < t - ?> 


( 20 ) 
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u 

x 


zfe (« - x) H(t • !>• 


( 21 ) 


Figure 24 shows the distribution of displacement, strain and velocity. 
When ct = a, all motion ceases and remains static, and the maximum 
distance the wave travels down the bar is not dependent on the 
characteristics of the bar but on the magnitude of the input force, 

T, and the magnitude of the friction force per unit length, F. 

Consider the static case where the force T is to the right and 
is balanced by F st a to the left, where F gt is the static function 
force. Since a, is defined as 



where F is the dynamic friction force, then 


F . (2T/ F) = T 
st 



This demonstrates that the static friction force is half the dynamic 
friction force; therefore, when the wave is traveling down the rod 
and reaches a , the friction force per unit length reduces from F to 
F/2, instantaneously. 

Also of interest in this problem is verifying the conservation 
of energy; i.e., the energy into the bar must be equal to the energy 
absorbed by the bar plus the energy loss due to friction. The energy 
relationship to be satisfied is the following: 
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Energy in = Kinetic Energy + Potential Energy + Friction Losses 
Tu = -| mu 2 + j Ku 2 + Fu 


Tu(0,t) 


PA 

2 



dt + f 


u^ dx + F 
x 



u dt. 


o 


o 


o 


With the appropriate substitutions from equations (19), (20) and 
(21), the following results; 

TFx q TFx 

4AE ^ X ° " = 4AE“ (x o " 2a ‘ ) ' 


This verifies the conservation of energy and partially verifies the 
correctness of the solution to this example. 

Of further interest is the conservation of momentum; i.e.. 



dt 



m du 


J ( T - Fct) dt = J pA ct du 
o o 


with the proper substitution from equation (19) the following 
results: 


Tt 



This verifies the conservation of momentum and also partially veri- 
fies the correctness of the solution to this example. 
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C. Square Wave Stress Impulse Loading to a Semi- Inf inite Rod 

This example, unlike the same example of the previous section 
without coulomb damping, may not be solved in the entirety with one 
solution, but this example must be subdivided for square waves of 
different pulse durations; i.e., t g aj c and t = aj 2c which are the 
only two square wave durations solvable by this technique. For one 
of these cases, a further subdivision of time is required. This will 
become apparent as these examples are analyzed. These subdivisions 
are required because of the inability to define a friction term which 
is all encompassing and, at the same time, is capable of being 
linearized. 

1. t g aj c 

This example is the case where a step stress is applied, and 
the wave travels along the rod until all motion ceases . Then at that 
time or some later time (depending on the value of t), the step stress 
at x = 0 will be released. This example will need to be further sub- 
divided: (1) For the time when the wave travels along the rod and 

comes to rest (t S t) and (2) for the time after the step stress is 
released (t > r). 

a. t S x 

This case has already been solved, and it is the step 
stress impulse loading of paragraph B of this section. 

b. t g x 

This case begins where the previous case ended; there- 
fore, the initial conditions of this case are the final conditions of 
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the previous case. By substituting into equations (19) and (20) 
t = a/c, the initial conditions are obtained; i.e., 

F 

u(x,0+) = (x - a ) 2 H(a - x) and u t (x,0+) = 0. 

Since it is known that the friction force in the static 
condition is F/2, and that the stress will be released causing a 
wave motion at that point with a velocity in the opposite direction 
to its previous motion, then the wave equation may be written as 
follows for this example: 


I? ' P SI* = • 2AE H(t ‘ f> + 2AE H(t) H(a ' x) * (22) 

Initially, the entire rod is under the influence of the static fric- 
tion force, but as time progresses and wave motion takes place, the 
friction force changes direction and doubles its value to meet the 
requirements of the wave motion. The boundary conditions at x = 0 
has changed from the previous case, because the stress is removed 
at x = 0; therefore, 

U x (0,t) = ° ° r U x (°’ s > = °* 

As x -> oo u(x,t) is bounded or U(x,s) is bounded. Taking the 
Laplace transform of equation (22), we obtain 
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Sfeil - U(x,s) + pr »(X,W) + ^ u t (x,CH-) 


2& e ’ lS + 2& Hto - X) - 


Substituting the initial conditions, using variation of parameters to 
obtain the particular solution and using the boundary conditions to 
resolve the constants, we obtain the following solution: 


-tse{ pr < t -! )2 + < 3 xc - 201 c)(t -t ) . 

+ (x 2 - 2ax + a 2 ) H(t) H(a - x)j- , 


H(t - f) 


(23) 


with the corresponding velocity and strain being 


u 

x 


" t • & <3ct • 2o ° H(t • 


(24) 


^ j(3x - 2a) H(t - |) - (2x - 2a) H(t) H(a - x)| . (25) 


Figure 25 shows the displacement, velocity and strain distributions. 
The velocity is similar to the velocity of the first part of this 
example except that all motion ceases at t = 2o/3c rather than 
t = al c. 

An energy check will be applied to this example to 
verify that energy is conserved. The Relationship to satisfy is 
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Original Potential Energy = Kinetic Energy + Potential Energy 

+ Friction Losses 

\ Ku 2 = mu 2 + Ku 2 + Fu 

a x o a x D 

J tu x (x,0)] 2 dx = J u 2 dx + ^ J u 2 dx + F J - (u - u Q ) dx. 

o o o o 

Substituting from equations (23), (24) and (25) using the proper 
forms, we obtain the following: 

F 2 ^ _ F 2 ^ 3 

24AE 24AE 

This verifies the conservation of energy and partially verifies the 
correctness of the solution to this example. The solution may be 
further verified by verifying the conservation of momentum; i.e.. 



o o 


With the appropriate substitution for equation (24), the following 
is obtained: 

3Fct 2 3Fct 2 

-Tt + , = -Tt + ■ , ■ . 

4 4 

This verifies that the solution did obey the law of conservation of 
momentum. 
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2. T = g/2c 

Consider the case of the semi-infinite rod with a square 
wave stress impulse loading vlrhere the pulse duration is t — O/J 2c. 

This case does not need to be subdivided into two time domains, 
because the friction force term can be defined for all time and be 
linear at the same time. It was observed when solving this case with 
the friction force acting only on the forward velocity wave (i.e., 
only on the region between ct - t § x § ct) that a negative velocity 
wave was generated immediately preceding the forward wave. There- 
fore, it was necessary to define a friction term which would resist 
the negative velocity rear wave as well as the positive velocity 
forward wave. This is done in the differential equation below: 


5 2 u 1 5 2 u F 

dx 2 " c 2 cit 2 AE 


H t 




The boundary conditions for a semi-infinite rod with a compressive 
square wave stress impulse loading are 


u(0,t) - - [H(t) - H(t - %)1. 


As x u (x, t) is bounded. Using these boundary conditions and 
zero initial conditions, we obtain the following solution: 



2x 


X 


a 


£ _L — _L — - — + ~ 

o: a 2 a 2 4 


( 26 ) 
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Figure 26 shows the distribution of displacement, velocity and strain. 
The front velocity wave and the rear velocity wave go to zero at the 
same time, t = a/c, and no motion exists after that time. Also, at 
t = aj c the strain has no discontinuity; indicating that no motion 
will be generated due to strain after that time. The static friction 
force per unit length that is required to maintain the strain is 



which is the same as that for the first part of the previous example 
with the exception that it is in two directions in this case with the 
change in direction at x = a/2. 

To verify the conservation of energy in this case, the energy 
balance must be subdivided into two time zones because of the dif- 
ficulty encountered in the required integration. Therefore, for 
t g t, the relationship to satisfy is 

x 0 x 0 x 0 

T u(0,t) = f uj dx + f J + F J u dx. 


o 


o 


o 
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With the proper substitutions, this results in 


2AE 


X “ 

T z x_ r x -\ 

2 - -2. 

°l 

1 

7 \ 

II 

a _ 

2AE |_ a J 


For t t, the relationship to satisfy is 



_c(t-T) ct 


c(t-x) ct 

- y> p 

T u(0,t) = Y 

j u 2 dx + / u| dx 

+ 2 

J a x dx *J U x dx . 


o C(t-T) 


o c(t-x) 


+ F 


_c(t-x) 


ct 


/ u dx + / " 


dx 


C(t-T) 


and with the correct substitutions, the following result is obtained 


3T g 0! _ 3T 2 o: 
8AE 8AE ’ 


which is a constant, as expected, because in this time zone no 
energy has been put into the rod. To verify the conservation of 
momentum, the same division in time is required. For t g t, 

t r 

J (T - Fct) dt = J pA-ct du 


Tt 


Fct 2 _ Tt _ Fct 2 

T ~ 2 ' 
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For t S t 


T t t U 

J (T - Fct) dt +J (-Fct) dt +J Fc(t - t) dt = J pA ct du 


u 

+ J' pA c(t - t) du. 


— - Fat + 


Fct* 


Ta 

= -r - Fct + 


With these verifications of conservation of energy and momentum, 
greater confidence in the solution is established. 

D. Step Stress Impulse Loading to the Left End 
of a Finite Rod with the Right End Free 

This example is a rod of length, £,, that has no supports at the 

right end and a compressive step stress is applied to the left end. 
As the wave traverses the tod, coulomb damping is present to retard 
the motion. The boundary conditions for this rod are 

- - aI H(t) or V°- s) - - aIs 

u(f,t) =0 or U(i5,s) = 0. 


The equation of motion for this case is simply equation (17) which 
is also the equation for the semi-infinite rod with a step stress 
input, because in both cases the velocity never changes direction. 
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With zero initial conditions and the above boundary conditions, the 
solution to equation (17) is 


-MZ[ 

n =0 


T(x - 2ni - x) - 7 (x Q - 2n.g - x ) 2 - ~ ( X Q 2n & " x ) 


4 v o 


oo — I 

H (t - && - f) [t(x q - 2 n^ + x) - & (x Q - 2 ni + x) 

n=l 
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+ ~ (x Q - 2 n^ - x) 


2nl + x ) + 

c c 


n=l 


| (x Q - 2nz - X ) 2 - ^ (x D - 2n$ - x) 


u i i- 2 n| _ x 
H l fc • “ cj\ ’ 


with the corresponding velocity and strain being 


n =0 


T + nFj e - 


Fx n 
o 




T - U 
1 2 


n=l 


H t - 


2 ni 


00 

-t) + X 

n=l 


T U 
2 


\--k{I[- T + T] 

n =0 

n=l 

To discover the limitations of these equations, we must re-examine 
the differential equation (17). The friction term for t > i/c 
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begins to apply beyond the length of the rod. While, in actuality, 
this poses no problem, it, nonetheless induces error into the mathe- 
matical analysis. To avoid this problem, the above equations will be 
limited to t g il c , and, in limiting the equations, they may be 
simplified to the following: 


u 


Ze{’ 


T(x c 


F 

*> + f 


(x £ 



u 


t 


u 

X 


AE 


AE 


i?-x -i 

T - -3»] H(t - f> 
-X +^H(t - |). 


(27) 

(28) 


(29) 


A new differential equation will be solved for t > c, viz. 


^ 2 u _ 1 d 2 u _ F 
dx 2 P St 2 AE 


v \ tt f *- \ 




Using the same boundary conditions as before but using end conditions 
of previous equations as the initial conditions to this analysis, 
i.e, , 


u(x, +0) = ^ 


T(l! - x) 


+ f (x 2 - 


Ji 2 ) and u fc (x, 40) 



we obtain the following solution: 
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\ / o x 2 2x 
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+ 1 l) + 4 V !? IP Jt 
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n=0 


^ 2n - 1 - - 2n - 1 - 5 


2 V^jg 


n=0 

with the corresponding velocity and strain being 


_ c r, w + v 

“t-sf ~ 2 + L 
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n=0 

H , , . iMiu . «)} 


n=0 

+ I[-¥]h(-^- 7 )}- 


(30) 


(31) 


(32) 


This solution is further limited in that T must be of sufficient 
magnitude that the velocity will never go to zero or negative at any 
point; therefore, T § FI!. The distribution of the velocity for 
T = YZ is obtained by substituting T = F<g into equations (28) and 
(31). For t S jg/c. 


u 


t 


Flc 

AE 



1 


(33) 
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For t § i/2, 


u -ns A 

t AE \2 


~ ^ [1/2] h (t - l 2 s±Ui + |\ 

n=0 ' 


00 



n=0 


H t 


(2n+l)j _ x^j_ 


(34) 


Equation (33) time, t l9 starts (i.e., t x = 0) when the step stress is 
applied. Equation (34) time, t 2 , starts (i.e., t 2 = 0) when the wave 
has traversed the length of the rod. Therefore, when t x = z/c, 
t 2 = 0. Figure 27 is the velocity for T = FZ. After the wave has 
traversed the rod once, the average velocity of the rod is constant 
(i.e., 5 t = 2AE ^ ’ therefore, the rod is traveling at a constant 
velocity while the wave motion is taking place. As t approaches 2Z/c, 
the velocity of the rod at x = 0 approaches zero, but at t < 2Z/c, the 
velocity is increased by the on-coming wave. Therefore, the velocity 
comes very close to zero, but never quite reaches it; this illustrates 
that T = FZ is the critical value. 

Now, consider a more typical value of T (i.e., T > FZ), e.g., 

T = 2FZ. Substituting into equations (27), (28), (29), (30), (31), 
and (32), results in the following, for t ^ i/c. 


u = ^ 
AE 


-2x 


2x 


,2 „ 


4 Z 2 4i 2 


H(t - f) 


_ FcJb 

U t AE 
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H(t - f>. 
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and for t § &lc. 


F4 2 J 3x o 7 _ 2x • 

“ aT1~ + 4 T + 
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x AE 
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n=0 


, , , (2n+l)£ , x 

[3/2] H 1 1 - c ' + ~ 


+ ^[-3/2] H (t - ( 2n t 1 I I - ^ 
n=0 


Figure 28 shows the displacement distribution; figure 29, the veloc- 
ity; and figure 30, the strain. The velocity and displacement are 
ever-increasing, while the strain has set up a cyclic pattern after 
the wave has initially traversed the length of the rod, the height 
of the discontinuity of the velocity and strain remains the same, 

|ip _ . The velocity at the boundaries increases each time by the 

discontinuity height upon reflection of the wave. 
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CHAPTER V 

Friction-Induced Wave Motion 

In this case, two rods are set together in such a manner that 
the wave motion in one rod induces wave motion in the other rod 
through friction (see figure 31). It is assumed that all cross- 
sectional planes remain plane, and the displacement is the same 
everywhere in the plane. The inducing rod will have a wave propa- 
gation velocity of c 1} and the induced rod will have a wave propa- 
gation velocity of c 2 . By using equation (17), the equation of 
motion becomes (for the induced rod): 


5 2 u 1 d 2 u F . 

S? ’ S? St* ‘ AE H<t 



Since these are semi-infinite rods and the induced rod has no strain 
at x = 0, then the boundary conditions for the induced rod are 

u(0,t) = 0. 

As x oo u(x,t) is bounded. Using these boundary conditions and 
setting the initial conditions to zero, the solution to the equa- 
tion of motion is 


c 2 c 2 

C 1 C 2 


U = 


2AE(c“ - c“ 
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) L C 1 c 2 
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with the following velocity and strain being 


u. = 


AE(cf 


2 2 

C 1 C 2 




) H(t - - (t - f-) H(t - f-) 
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l gi c i 

AE(c 2 _ c 2) 


(t 


H(t - 7") + (t - f-) H(t - 
c 2 



To gain an understanding of this example, three cases will be 
examined: subsonic, “ 2 = sonic, = 1 and supe sonic, — 1 ~ - 2. 
Figure 32 shows the distribution of displacement, velocity and strain 
for these three cases. Observe the similarity in the wave shapes of 
the subsonic and supersonic cases, and the difference in the sonic. 

In the subsonic case, the wave has propagated twice the distance down 
the rod as does the sonic and supersonic case because, in the sub- 
sonic case, the limiting factor is the inducing rods velocity, c x ; 
and in the sonic and supersonic case, the limiting factor is the 
induced wave velocity, c 2 . The displacement is greater for the cases 
in this order: subsonic, sonic and supersonic, because the friction 
force is acting over a greater portion of the rod for a longer period 
of time in that order. The strain is greater in the same order for 
the same reason. 
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CHAPTER VI 

Fourier Transform Technique 

In solving examples similar to the ones in this treatise, the 
Fourier transform is taken of the axial coordinate variable, x, 
rather than the time variable, t, as was the case in Laplace trans- 
forms. Instead of having initial conditions as associated with the 
Laplace transforms, the strain boundary condition at x = 0 and for 
all time is required. The initial conditions are needed to resolve 
the constants of integration in the same way that the boundary condi- 
tions were used in Laplace transform technique. Therefore, in solv- 
ing these problems, three quantities must be known: (1) u x (0,t), 

(2) u(x,0), and (3) u t (x ? 0). It now becomes apparent that the 
Fourier transform is limited to solutions of examples with semi- 
infinite rods, because there is no way to stipulate the boundary 
condition at x > 0 .. 

A. Semi-Infinite Rod with a Step Stress Impulse Loading 

For a compressive step stress impulse loading, the boundary 
condition at x = 0 is a constant; i.e., 

u (0,t) = - -t H(t) . 
x' ’ ' AE 


Using the wave equation (2) and taking the Fourier transform yields 
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-p2u(p,t> -yo.t) 

+ c*p* U (p>t) - giH(t) 
which has a solution of 

U(p,t) = B!e 1Cpt + B 2 e _:LCpt + H(t). 

Since the initial conditions are zero; i.e., 


u(x, +0) = 0 and u fc (x, +0) = 0, 


U(p,t) 



cos 

n 



( 35 ) 


Because it is very difficult to obtain the solution by integration 
of the inverse Fourier transform integral and because insufficient 
Fourier transform tables are available, the solution will be proved 
by taking the Fourier transform of the known correct solution. The 
correct solution from chapter III is 


u = H(t-.f). 


Substituting this equation into the transform integral yields 
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OO 

U(p,t) = ~tkf (ct - x) H(t - cos px dx 


ct 


U(p,t) = 


AE 


ct 

ct cos px dx - J' x cos px dxj- 


Upon completing the indicated operation and some simplification of 

/ 

terms, we obtain 


U(p 


T fl - cos cptl 
•*> * AE { ’ 


which is identical to equation (35). 

B. Semi-Infinite Rod with a Square Wave Impulse Loading, t = aj 2c 

This example must be subdivided into two time domains; 
t I o/' 2 c and t ^ aj 2 c. 

For t % aj 2c, the boundary condition is 


= - AE 


and the initial conditions are zero. Using these conditions and 
equation (17) results in the following: 


U(p,t) = 77 ft (1 - cos cpt) + | cos cpt - S 1 p 3 Cpt ^ )| • 


AE p 


( 36 ) 
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The correct solution is equation (26) when it is limited to t S OtJ 2c, 


u 



Tx 


F c 2 t 2 
4 



Taking the Fourier transform of this equation, gives the following 
results : 


U(p,t) “ n {f (1 “ COS Cpt) + 2 


ct 1 

—n cos cpt H — 3 sin cpt 

p P 


]} 


which is identical to the solution obtained in equation (36). 

For t g a/ 2c, the boundary condition is 

u(0, t) = 0 

and the initial conditions are the end conditions of the above case 
i.e. , 


The differential equation is based on a friction force that already 
exists from cd 2 to 0 and moves across the rod with the wave, and 
a negative friction force to resist the negative rear velocity wave 


which follows the forward wave, i.e., 
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d 2 u 1 d 2 u _ F 
c)x 2 ^ ^t 2 AE 


H (ct + | - x) - 2H(ct - x) 


Using the aforementioned boundary and initial conditions, taking 
the Fourier transform, and doing much simplification of terms, we 
obtain the following solution: 



Again, equation (26) is the correct solution when t is set equal to 

t + cc/ 2c, because this solution began when equation (26) was at 

CC 

ajlc. Therefore, substituting t = t + — into equation (26) yields 


rr 


u = — i'i T f ct 

AE \l \ 


a \ Fc^r^ 

2 ■ X J ■ ~T~ 


Fact _ vor , Fx^ 

4 16 4 


H 



T(ct - x) 


Fc 2 t 2 Fx 2 ' 
2 2 


H(ct - x) 



Taking the Fourier transform yields 


U(p,t) 



cos cpt - cos 


(cpt + 


Fct 

IP 


2 cos cpt 


- cos I cpt + 


232. 

2 


')] + 2? ; 


2 sin cpt - sin ( cpt + 


»)]} 


which is identical to equation (37). 



CHAPTER VII 


Summary and Conclusions 

The applicability of Laplace transforms to aid in solving wave 
motion differential equations, excluding damping, was demonstrated 
in that many types of boundary conditions were readily solved. For 
this wave motion without damping, many interesting characteristics 
were exhibited, especially at the boundaries of the rod. Of special 
interest are figures 7 , 11 , 12 , 17 , 18 and 20 . The complexity of 
the wave shapes at the boundaries exhibits the details of the wave 
motion characteristics, and from these details it becomes apparent 
that difficulty will be encountered in defining a friction force 
which will always retard the motion and at the same time be linear. 

In this treatise the interesting characteristics of wave motion 
with coulomb damping occur along the length of the rod and not at 
the boundaries. Perhaps the boundaries would have proved to be more 
interesting if more complex solutions had been obtained. In all of 
these solutions, it was found that the velocity was a function of 
time only and not the axial distance except for the sudden change at 
the point of discontinuity. The strain was a function of the axial 
distance only except for the shifting of the point of discontinuity 
with time. The displacement was a function of time and axial 


distance . 
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For the semi- inf inite rod with the step stress impulse loading 
and coulomb damping, the maximum distance that the wave would pro- 
pagate along the rod ( a = 2T/F) was not a function of the charac- 
teristics of the rod, but was dependent only on the magnitude of the 
input force and the friction force per unit length. The static 
friction force per unit length that existed after all motion ceased 
was one-half that of the dynamic. When the step stress was removed, 
the maximum distance the wave propagated was 2a/ 3. The static 
friction force between 0 and 2o/3 was one-fourth the dynamic fric- 
tion force, and between 2a/ 3 and a, it was one-half. For the same 
rod with a square wave impulse loading of pulse duration x = aj 2c, 
all motion ceased at a, and the static friction force was one-half 
the dynamic friction force. 

For the finite rod with a step stress impulse loading and 
coulomb damping, it was found that for T = F£ y Lhe average velocity 
of the rod was a constant (Fc y g/2AE) after the wave has traversed the 
rod once while wave motion continued. For T > F Jb 9 the velocity con- 
tinues to increase while the strain sets up a cyclic behavior. 

Once the wave has traversed the rod, the discontinuity height of 

i F£ i 

the velocity and strain is a constant, |T - -~| , and the velocity 
at the boundaries is increased by the amount of the discontinuity 
height (T j-) upon reflection of the wave. 
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VELOCITY AND STRAIN DISTRIBUTIONS 
TE ROD WITH A STEP STRESS 
>ULSE LOADING 
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STRAIN 



FIG. 3. DISPLACEMENT, VELOCITY AND STRAIN DISTRIBUTIONS 
OF A SEMI-INFINITE ROD WITH A SQUARE WAVE STRESS 
IMPULSE LOADING FOR tMj. 
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FIG. 4. DISPLACEMENT, VELOCITY AND STRESS DISTRIBUTIONS 
OF A SEMI-INFINITE ROD WITH A STEP VELOCITY 

IMPULSE LOADING 
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FIG. 5. MATHEMATICAL MODELS REPRESENTATION OF 
THE DISPLACEMENT, VELOCITY AND STRAIN OF 
A ROD WITH A STEP STRESS IMPULSE LOADING ON 
THE LEFT END AND FREE ON THE RIGHT END AT l»0 






FI6. 6. DISPLACEMENT, VELOCITY AND STRAIN DISTRIBUTIONS OF 
A ROD WITH A STEP STRESS IMPULSE L0ADIN6 ON 
THE LEFT END AND FREE ON THE RIGHT END 
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FIGURE 7 (b) STRAIN DISTRIBUTION OF A ROD WITH A SQUARE WAVE STRESS IMPULSE L0ADIN6 

ON THE LEFT END AND FREE ON THE RI6HT END, r * //2c 
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FIG. 8. 
IMPULSE 


DISPLACEMENT DISTRIBUTIONS OF A ROD WITH A SQUARE WAVE STRESS 
LOADING ON THE LEFT END AND FREE ON THE RIGHT END, r * //2c 
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FIG. 10 VELOCITY 8 STRAIN DISTRIBUTIONS OF A ROD WITH AN ATTACHED STEP 
VELOCITY IMPULSE LOADING ON THE LEFT END 8 FREE ON THE RIGHT END 



FIG. 11. VELOCITY AND STRAIN DISTRIBUTIONS FOR 
A ROD WITH A SQUARE WAVE VELOCITY IMPULSE 
LOADING {t * ^r) ON THE LEFT END AND FREE ON THE RIGHT END 
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THE SAME THING HAPPENS AT X * 0 
AND FROM THEN ON, AT EACH END 


FIG. 13. VELOCITY DISTRIBUTIONS OF A ROD WITH A TRIANGULAR STRESS WAVE 
IMPULSE LOADING ON THE LEFT END AND FREE ON THE RIGHT END 
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FIG. 14. STRAIN DISTRIBUTIONS OF A ROD WITH A TRIANGULAR STRESS WAVE 
IMPULSE LOADING ON THE LEFT END AND FREE ON THE RIGHT END 





FIG. 15. DISPLACEMENT DISTRIBUTIONS OF A ROD WITH A TRIANGULAR WAVE STRESS 
WAVE IMPULSE LOADING ON THE LEFT END AND FREE ON THE RIGHT END 
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FI6. 18. DISPLACEMENT DISTRIBUTIONS OF A ROD WITH A SQUARE WAVE STRESS 
IMPULSE LOADING ON THE LEFT END AND FIXED ON THE RIGHT END 
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FIG. 19. TRIANGULAR IMPULSE LOADING WAVE 
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FIG. 20. VELOCITY AND STRAIN DISTRIBUTIONS OF A ROD WITH A TRIANGULAR 
WAVE STRESS IMPULSE LOADING ON THE LEFT END AND FIXED ON THE RIGHT END 
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FIG. 21. DISPLACEMENT DISTRIBUTIONS OF A ROD 
WITH A TRIANGULAR WAVE STRESS IMPULSE LOADING 
ON THE LEFT END AND FIXED ON THE RIGHT END 


>o|x >*|x 



FI6. 22. WAVE MOTION WITH COULOMB DAMPING 

FREE BODY DIAGRAM 



FIG. 23. WAVE MOTION WITH COULOMB DAMPING 
AND STEP STRESS IMPULSE LOADING 
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FIG. 24. DISPLACEMENT, VELOCITY AND STRAIN DISTRIBUTIONS 

FOR A SEMI-INFINITE ROD WITH COULOMB DAMPING 

AND A SQUARE WAVE STRESS IMPULSE LOADING FOR t > « AND t < r 

c 













FIG. 26. DISPLACEMENT, VELOCITY AND STRAIN DISTRIBUTIONS 
FOR A SEMI-INFINITE ROD WITH COULOMB DAMPING 
AND ASQUARE WAVE STRESS IMPULSE LOADIN6 OF T/A FOR r*£ 






FIG. 28. DISPLACEMENT DISTRIBUTION 
FOR A FINITE ROD WITH COULOMB DAMPING 
AND A STEP STRESS IMPULSE 

LOADING OF 4 - FOR T*2F 1 








FIG. 29. VELOCITY DISTRIBUTION FORA FINITE ROD 


WITH COULOMB DAMPING AND A STEP STRESS 
IMPULSE LOADING OF t FOR T = 2 F i 
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FIG. 30. STRAIN DISTRIBUTION 
WITH COULOMB DAMPING AN 
IMPULSE L0ADIN6 OF £ 
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FI6. 32. DISPLACEMENT VELOCITY AND STRAIN DISTRIBUTIONS FOR FRICTION - INDUCED WAVE MOTION 
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